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We investigate flow separation in the air flow over transverse sand dunes. CFD simulations of the air flow 
over differently shaped dunes are performed. The length of the recirculation region after the brink of the dune 
is found to depend strongly on the shape of the dune. We find that the nondimensionalised separation length 
depends linearly on the slope of the dune at the brink within a large interval. A phenomenological expression 
for the separation length is given. 
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1 INTRODUCTION 

Dunes are naturally occurring, beautifully shaped 
sand deposits. Since the middle of the previous cen- 
tury, they have attracted the attention of scientists who 
have been seeking to model them and understand the 
processes leading to their formation. From the point 
of view of the physicist, sand dunes constitute a vari- 
able boundary problem: The air flow is determined by 
the shape of the dune and in turn influences the dune 
shape by transporting sand grains. Therefore the air 
flow over dunes is of great importance for understand- 
ing dune formation and evolution. 

Since the start of scientific interest in 
dunes, there has been some work on this 
topic, both theoretical (INelson and Smith~1 989: 
Parsons, Walker, and Wiggs 200 4|) and ex peri- 
mental ( |Engel 1981| ISweet and Kocurek 1 990). 
However, due to the difficult nature of the prob- 
lem, these works have only tackled part of the 
problem. In several publications, transverse 
dunes have been modeled as triangular struc- 



tures (Engel 198 11 Parsons, Walker, and Wiggs 2004 
|Parsons, Wiggs, Walker, Ferguson, and Farvey 2 004 ) 
Field measurements of air flow over 
dunes tend to lack measureme nts of the 

dune profile (ISweet and Kocurekl 990: 

Frank and Kocurek 19961) . 

A recent field measure- 



the foot of the following dune shape is located 
at or downwind of the flow reattachment point. 
Therefore the distance of closely spaced dunes is a 
measure of the length of the recirculation region. 



In dParteli, Schwammle, Herrmann, Monteir o, and Maia~ 2004 
the separation length after different dunes was deter- 
mined in this way. 

In this work we will present results for widely 
spaced or isolated transverse dunes. This is to some 
extent an idealisation. However, we think this is valid 
and useful: We want to concentrate on the effect of the 
dune shape, other things being equal, and the presence 
and shape of neighbouring dunes would constitute ad- 
ditional parameters. 

2 METHOD 

Our simulations were performed with the computa- 
tional fluid dynamics software FLUENT. The wind 
flow over dunes is fully turbulent. We simulated the 
Reynolds-averaged Navier-Stokes equations using the 
k-e model for closing the equations. 

The simulations were two-dimensional, implying 
a wind direction perpendicular to the dunes. The 
cross sections of the dune shapes were constructed 
from two circle segments, a concave one mod- 
eling the foot of the dune and a convex one for 
the crest. We verified that this is realistic enough 
to describe the transverse dune shapes measured in 



ment (Parteli, Schwammle, Herrmann, Monteiro, and 
suggests that the shape of transverse dunes has sig- 
nificant influence on the length of the recirculation 
region. Since the sand transport in the recircu- 
lation region in the lee of a dune is negligible, 



Mai(PM©I) Schwammle, Herrmann, Monteiro, and Maia 2004) 
well. To obtain different shapes, the position of the 
slip face was varied from the start to the end of the 
convex part, see Figure [lj Note that this has the 
consequence that not all the dunes have the same 
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Figure 1: The seven different dune shapes investi- 
gated. The scale displays the brink position. The 
crest height of the dunes with negative brink position 
equals the brink height, the height of those with posi- 
tive brink position is 3 metres. 
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Figure 2: The geometric variables characterising the 
dune shapes. The brink angle is positive for dunes 
with a sharp brink and negative for round dunes as 
the one shown in the figure. 

height. The respective heights and other geometrical 
data are given in Table [lj 

The velocity profile at the influx boundary of the 
simulation region was set to the logarithmic pro- 
file which forms in flow over a plane in neutral 
atmospheric conditions. The shear velocity was 
chosen to be 0.4 m/s, the roughness length 25/im. 
The size of the roughness elements on the ground, 
the grain size, was set to 250yum. The roughness 
length is generally considered to be of the order 
of the grain size or up to a factor of 30 smaller 
( |Wright, Schaffner, andMaa 1997 ). 

The region around the dune in which the flow was 
simulated was chosen large enough so that the bound- 
aries did not influence the results. This was verified 
by choosing larger simulations areas for some dune 
shapes and comparing the results. From the brink po- 
sition 0, the simulation region extends 45 m to the left 
and 70 m to the right (see Figure[3]). The height of the 
simulated region was chosen to be 30 m for all dunes 
except the one with brink position -15 m, where 20 m 
was found to be sufficient. 



of the simulation grid. Therefore we performed the 
simulation of the flow over each dune with three dif- 
ferent grids and interpolated the separation lengths to 
the continuum. The average grid spacings were 10, 7 
and 5 cm, respectively. 

3 RESULTS 

Table [T] shows the results for all dune shapes. The 
length of flow separation, our quantity of interest, was 
measured from the slip face brink, where the flow sep- 
arates, to the flow reattachment point (see Figure |2j). 
The errors were estimated to be one grid spacing for 
the determination of the flow reattachment point. The 
separation lengths determined for the different grids 
and their errors were interpolated to the continuum 
with the standard linear regression formulas. To the 
statistical error we added quadratically a systematic 
error of 0.5 metres. The systematic error accounts 
for biases which may be inherent in the turbulence 
model and parameter settings used. We expect it to be 
strongly correlated. 

To nondimensionalise the separation length /, it 
was divided by the height of the slip face. This di- 
mensionless quantity is universal, that is it does not 
depend on the absolute height of the dune. That is 
because the flow is already fully turbulent for small 
dunes. We find that //A is larger for dunes with a 
sharp brink than for rounded dunes. It depends lin- 
early on the brink position d respectively the angle 
of the dune shape at the brink, a. As can be seen in 
Figure El the linear relation extends up to an absolute 
angle of 7.5°. This amounts to a variation of the dune 
length between 6.7 and 13.3 times the height. Fitting 
the relation 



Z(a)/A(a) = A-a + B. 



(1) 



we obtain A = 0.208 /° and B = 5.85. The points with 
brink position ±15 were ignored for this fit since they 
deviate from the linear law. 

Since the brink angle a and the brink position d are 
related by the geometry of the dune, the fit ® can be 
reformulated to give the separation length in terms of 
the brink position: 



30 m 



115 m 



Figure 3: The simulated region around the dune. 

The length of flow separation, our quantity of in- 
terest, was found to depend slightly on the spacing 



l(a(d)) = (A ■ a(d) + B) A(a(d)) 



d 



-A arcsin — + B 
R 



1 



(2) 



• ^max - d tan arcsin -^J j 

This equation contains the height of the round dunes, 
-f^max = 3 m, and the radius of the circle segments used 
to model the shape, R = 75.75 m. Remarkably, this 
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Brink pos. 
d [m] 


Height 
H[m] 


Brink 
height A [m] 


Angle at 
brink a [°] 


Separation 
length I [m] 


Error 
estimate [m] 


l/A 


-15 


1.5 


1.5 


11.4 


11.74 


±0.026 ±0.5 


7.83 


-10 


2.337 


2.337 


7.6 


17.33 


±0.048 ±0.5 


7.41 


-5 


2.835 


2.835 


3.78 


18.88 


±0.042 ±0.5 


6.66 





3 


3 





17.30 


±0.065 ±0.5 


5.77 


5 


3 


2.835 


-3.78 


14.44 


±0.050 ±0.5 


5.09 


10 


3 


2.337 


-7.6 


9.97 


±0.041 ±0.5 


4.27 


15 


3 


1.5 


-11.4 


4.96 


±0.021 ±0.5 


3.31 



Table 1 : Geometry of the simulated dunes and results for length of flow separation. See Figure[2]for a definition 
of the geometric variables. The brink angle is defined to be positive if the upwind slope is positive at the 
brink. The first error in the separation length is the statistical error in the determination of the length from the 
simulation data, the second error is the systematic error of the simulation (see text). 
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Figure 4: Dependence of the nondimensionalised sep- 
aration length on the angle a. The relationship is re- 
markably linear. The rightmost value of a belongs to 
the dune with the sharpest brink. 
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Figure 5: Dependence of the flow separation length 
on the brink position. The expression derived from 
the linear angle dependence displayed in Figure|4]pro- 
vides a much better fit than a parabola. 



expression fits significantly better than a parabola (see 
Figure |3) even though it has a smaller number of pa- 
rameters. It has a x 2 = 0.36 compared to x 2 = 4.3 for 
the parabola fit. 

Equation |2] can be approximated by using 
arcsinx ~ x and tana; ~ x. In this way one obtains a 
polynomial expression for the dependence of the sep- 
aration length on the brink position: 



1(d) 




^d 3 - — d 2 - 
2 R 2 2R R 



Cd s + Dd 2 + Ed + F 



d + BH n 



(3) 



The coefficients are: C = 0.00103 m" 2 , D = 
0.0386 m" 1 , E = 0.472, and F = 17.55 m. The val- 
ues of this expression are indistinguishable from © 
in the range of d which we investigated. Therefore 
we do not plot it separately in Figure 

4 DISCUSSION 

Comparison of our results to other work is hampered 
by the fact that the dependence of flow separation 
on the dune shape has not been investigated before. 
Therefore, the following comparisons are to be un- 
derstood as consistency checks. 

A recent review of air flow over transverse dunes 



(Walker and Nickling 2002) cites values of 4-10 for 
l/A. Our results also lie within that range (see 
Figure |4|). Engel (Engell981) finds values for 
the nondimensionalised separation length between 
4 and a little over 6, depending on the rough- 
ness and the aspect ratio of triangular dunes. In 
( Parsons, Walker, and Wiggs 2004| ) a wide range of 
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between 3 and 15 is given for the same quantity. The 
values for an aspect ratio of 0.1, which applies to 
our dune with a = 0, are 5.67 and 8.13, depending 
on the height. This compares well with our value of 
5.76. The discrepancy can be explained by the differ- 
ent shape, in particular the fact that our dune shape for 
a = has a hor izontal tangent at the brink, whereas 
the dunes in ( |Parsons, Walker, and Wiggs 2004| ) are 
triangular. 

Unlike our simulations, which treated an 
isolated dune shape, the field measurements 



(Parteli, Schwammle, Herrmann, Monteiro, and Maia 
were performed in a closely spaced dune field. The 
authors find that the distance between the brink 
of each dune and the foot of the following one is 
typically four times the height or below. Under 
the assumption that the dune field is station- 
ary, this distance is a measure of the separation 
length. This is at the lower end of the separa- 
tion lengths we obtain, for very rounded dunes. 
The dunes with the smallest separation length in 
( Part eli, Schw ammle , Herrmann, Monteiro, and M aia 
are indeed very round. The small values for other 
dunes can be put down to the different situation of 
closely spaced transverse dunes. 

5 CONCLUSIONS 

We have determined the length of flow separation 
in the lee of isolated transverse dunes of different 
shapes. The separation length nondimensionalised by 
division by the slip face height was found to be larger 
for dunes with a sharper brink. For a wide range of 
dune shapes, this growth is well described by the lin- 
ear relationship ([T]). 

The dependence of the absolute separation length 
on the position of the slip face is described by the ex- 
pression © more accurately than by a parabola. The 
maximal separation length does not occur for dune 
shapes with a horizontal tangent at the brink, but for 
shapes with a somewhat sharper brink. 

These results were obtained for isolated dunes up 
to 3 metres in height. It would be interesting to see 
how the flow separation length is influenced by the 
presence of other dunes nearby, as in a closely spaced 
dune field. This topic will be treated in a future pub- 
lication. 
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